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We study the interaction of a liquid drop with an elastic beam in the case where bending effects dominate. 
We use a variational approach to derive equilibrium equations for the system in the presence of gravity 
and in the presence or absence of contact Une pinning. We show that the derived equilibrium equations 
for the beam subsystem reveal the external forces applied on the beam by the liquid and vapor phases. 
Among these, the force applied at the triple hne (the curve where the three phases meet) is found to he 
along the liquid-vapor interface. 
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1. Introduction 

The present trend toward miniaturization of engineering systems and machines is giving surface 
effects the leading role: in a system of size L the respective scalings of volume (L^) and surface 
(L^) forces are such that the latter outrange the former as soon as L is small enough. Among 
other surface effects, surface tension is now widely used at small scales, for example to self- 
assemble microsystems (Syms et al, 2003). The concept of force is not easy to explain to 
recalcitrant students: Has anyone already seen a force? How to be sure of the direction of an 
applied force? Capillary forces are no exception and conceptual questions about it recurrently 
emerge (Marchand et al, 2011). As these capillary forces are now used to bend small elastic 
structures (Roman and Bico, 2010) and as it has been recently proposed that these forces might 
not act as previously thought (Marchand et al, 2012), we here study the interaction of a liquid 
drop with a flexible beam from an energy point of view: we derive equilibrium equations of 
the system from a variational approach that is merely built on the classical hypothesis of the 
presence of surface energies arising at interfaces between the three phases: solid, liquid, and 
vapor. 

In section 2 we recall that the Young-Dupre relation for the contact angle of a drop lying on 
a substrate can be derived from a variational approach where the concept of force is not invoked, 
as first realized by Gauss (1830, 1877). In section 3 we consider the case where the substrate is a 
flexible beam and we add gravity (for the beam and the liquid) in Section 4. We recall in Section 
5 that contact line pinning can also be considered from an energy point of view, and we finally 
illustrate our model by computing the behavior of a drop-beam system as the drop evaporates. 
Conclusion follows in Section 6. 
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2. Liquid drop on rigid substrate, no gravity 

We consider the equilibrium of a liquid drop of given volume sitting on a rigid substrate of length 
L, see Fig. 1. If the drop is small enough, gravity and the hydrostatic part of the pressure can 
be neglected and consequently the liquid-air interface is circular. For simplicity we adopt the 
two dimensional framework introduced in Rivetti and Neukirch (2012) where the liquid- vapor 
interface is a cylindrical arc; We call r its radius, w its height, and 2/3 its opening angle. The 
liquid vapor interface then comprises (z) a cylindrical surface of area 2/3rw and (ii) two planar 
caps, each of area A = r^(/3 — sin (3 cos /?). The wetting angle is equal to /3 and the wetted length 
of the beam is noted 2D. To each of the three different interfaces, liquid-solid, liquid- vapor, and 
solid-vapor we associate an energy per area: , 7^^, and 7^^ respectively. The energy of the 
system is then given by the sum: 

E{(3, r, D) = 2{w^r + A)jev + 2wD-fes + 2w{2L - D)-fsv (2.1) 

To minimize E under the constraints of (i) fixed volume V = Aw = wr'^{j3 — sin /3 cos /3) and 




Figure 1. A two-dimensional model of a liquid drop at rest on a rigid substrate. In the absence of gravity the liquid- 
vapor interface is a circular arc of radius r. The contact angle /3 is set by the balance of the interfaces energies, see 
Eq. 2.5, and the wetted length 2D depends on the drop volume V. 

{a) geometric relation D = r sin /?, we then introduce the Lagrangian: 

C = E- r]w{2D - 2r sin /3) - [r'^{l3 - sin /3 cos (3)] (2.2) 

where rj and n are Lagrange multipliers. We note 7 = 7^ and A7 = 7^^ — Jsv- Equilibrium 
equations are found by imposing that derivatives of C with regard to the three variables (3,r,D 
vanish: 

— =0 = 2wrj + 2r]wr cos f3 — fiwr'^{l — cos 2/3) (2.3a) 
dC 

— — =0 = 2w(3'y + 2'nw sin (3 — 2u,wr{(3 — sin (3 cos (3) (2.3b) 
or 

dC 

=0 = 2u;A7 - 2r]w (2.3c) 

oD 
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Combining (2.3a) cos/3 + (2.3b) r sin/3 and (2.3a) sin/3 - (2.3b) r cos/3 yields 'j = fir and 
7 cos /3 + 77 = 0, and using (2.3c) gives r] = A7. We finally arrive at: 

A7 + 7 cos /3 =0 (2.4) 

M=- (2.5) 

r 

The first equation is the well-known Young-Dupre relation giving the contact angle and can be 
interpreted as a force balance of the triple point in the horizontal direction. The second equation 
gives the Laplace pressure inside the liquid drop. In the vertical direction force balance is also 
achieved: the vertical forces acting on the rigid substrate are the distributed Laplace pressure jj, 
and the surface tension "f£y, with the total downward force being 2Dwfi and the total upward 
force being 27^j, sin /3. Using (2.5) and D = r sin /3, we see that these two forces equilibrate. We 
shall see in the next section that when the substrate is a thin elastic strip, these forces induce 
flexural deformations. 



3. Liquid drop on a flexible beam, no gravity 

We now consider the case of a liquid drop sitting on an elastic strip, see Fig. 2, and we look 
for equilibrium equations governing the bending of the elastic strip by capillary forces. We 
still work under the hypothesis where gravity and the hydrostatic part of the pressure can be 
neglected, yielding a circular liquid-air interface. In addition to the sum of the three interface 
energies: 

= 2{w/3r + A)7to + 2wDjes + 2w{2L - D)^sv (3.1) 

we consider the bending energy of the elastic strip. We use the arc-length s along the strip 
to parametrize its current position r{s) = {x{s),y{s)). The unit tangent, t{s) = dr/ds, makes 
an angle 9{s) with the horizontal axis: t = {cos 9 {s),sui 9 {s)). The bending energy density is 
proportional to the square of the curvature 9'{s)\ 




Figure 2. An elastic strip bent by capillary forces. In the absence of gravity the liquid- vapor interface is a circular arc 
and the strip outside the interaction region is straight. The system is invariant is the z direction, with width w. 
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[B'{s)fds 



(3.2) 



where YI is the bending rigidity of the strip {Y is Young's modulus of the beam material, 
and I = h^w /12 is the second moment of area of the section of the beam). We minimize E = 
+ E^ under the following four constraints. First, the liquid volume V is fixed. It is given by 
V = wA, where A is the area in between the liquid-air interface and the liquid-solid interface: 



A = r^ I /3- ^sin2/3 ) + 2x{D)y{D) - 



•x(D) 
x{-D) 



ydx 



(3.3) 



Second, we have the geometric constraint: 



x{D) = rsmP 



(3.4) 



that is due to the intersection of the circular liquid- vapor interface and the beam at s = D. As 
the variables x{s), y{s), and 6{s) all appear in the present formulation, we have to consider the 
continuous constraints relating them. These are our third and fourth constraints: 



x'{s) = cos9{s) , y'{s) = sin6{s) 



(3.5) 



These continuous constraints necessitate the use of varying Lagrange multipUers u{s) and A(s). 
We therefore introduce the Lagrangian: 



£ =E^ + Ej — ijlw 

— r]w [x{D) — r sin /3] + 



rD 



- -sm2pj +2x{D)y{D)- ^ yx'ds 



■L rL 

z/(s) [x' - cos 9] ds + A(s) [y' - sin 9] ds (3.6) 

-L J-L 



As we will only treat symmetric solutions, we focus on the positive s interval: s G [0; L] with the 
following boundary conditions: 



x{0) = 0, y(0) = 0, 9{Q) = 



(3.7) 



We also remark that in this elastica model, the varying Lagrange multiplier and A(s) will 
be found to be the internal force components, see Appendix 1 . As the external force coming for 
the meniscus will make the internal force discontinuous as s passes through s = D,v/e split the 
two last integrals in the Lagrangian (3.6) and write: 



sG[0;D) : i^{s) = i^iis) , \{s) = Xi{s) 
s G {D; L] : i^{s) = i^e(s) , A(s) = Ae(s) 



(3.8) 
(3.9) 
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C{x, y, 9, p, r, D) = -YI [e'{s)] ^ ds + w/3r-i + wDA-f - ■qw 



+ 



+ 



•D 




f-D 



x'ds — r sin /3 



Uo 



1 

2' 





rD 






+ 


x'ds X 


y'ds - 


yx' ds 







J 














i'i{s) \_x' — cos^] ds + fe(s) \_x' — cosO\ ds 



[y' - sin 6*] ds + Ae(s) [y' - sin 6'] ds 

JD 



(3.10) 



where 7 = 7^ and A7 = 7^5 - 7^ 



Fjrj? variation 

The energy £' and the Lagrangian C are functions of the variables x, y, 9, P, r, D. We note 
X = {x,y,9, /3, r, D) and we consider the conditions for the state to minimize the energy E. 
Calculus of variation shows that a necessary condition is: 



C'{X,){X)=-C{X, + eX) 



(3.11) 



e=0 



where X = (,t, y, ^, /3, r, i^). Moreover boundary conditions (3.7) implies that x(0) = 0, ?/(0) 
0, 6'(0) = 0. Noting that: 



[■A+eA rA 

f{x)dx = f{x)dx + €Af{A) + 0{e^) 
Jo 



(3.12) 
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we evaluate the first variation (3.11) to be: 



C'{X,){X) = YI 



e' e'ds + w/3f-f + iw/3r7 + wDAj 







— fiw 



rr [ 0-^ sm2p\ + ^ (1 - cos 2^) 







rD 




rD 






x'ds X 


y'ds + 


x'ds X 
















D 



y'ds + Dx'{D)y{D) + Dx{D)y'{D) 



— ryw 



■D 




yx' ds + Dy{D)x'{D) 



x'ds + Dx'{D) — rsm/3 — rP cos (3 



+ 



+ 



rL 







Ui{s) [x' + e sin 6'] ds + Ue{s) [x' + 6 sin 9] ds 



D 



Xi{s) [y' - 9 cos e] ds + Ae(s) [y' - 9 cos 9] ds 

JD 



(3.13) 



where we have used (3.5) at s = D to eliminate sorne terms related to the last four integrals. We 
require this expression to vanish for all x{s), y{s), 9{s), P, f, and D. For (3.13) to vanish for all 
/3, we must have, as before: 



rw-f — l^w— (1 — cos 2/3) + rjwr cos /3 = 



(3.14) 

For (3. 13) to vanish for all f, we must have, as before: 

Pwj — iiwr {(3 — sin/3 cos/3) + rjw sin (3 = (3.15) 

Combining these last two equations we obtain: 

^r = 7 and 7cos/3 + ?7 = (3.16) 

where /j, is identified to the Laplace pressure. For (3.13) to vanish for all D, we must have: 

wA-f - iJ.wx{D)y'{D) -r]wx'{D) = (3.17) 

Using (3.4), (3.5), and (3.16) we obtain: 

A7 + 7cos[/3 + ^(i:>)]=0 (3.18) 

This is the Young-Dupre relation for the wetting angle /? + 9{D) between the beam and the 
Uquid-air meniscus. Requiring (3.13) to vanish for all 9 yields, after integration by parts: 



YI[9'9]^+ [-YI9" + i^i sin 9 - Xi cos 9] 9ds 
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+YI[9'9]^+ [-YI9" + i^esm 9- Xe cos 9] 9ds = (3.19) 

JD 
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which implies that the curvature 6'{s) is continuous as s goes through D and that it vanishes at 
the s = L extremity. Moreover we obtain the moment equilibrium equations: 



Yie" =Ui sin e - Xi COS for s e [0; D) 
Yie" =Ue sin e-\e cos e iorse{D- L] 



(3.20) 
(3.21) 



where we see that the continuous Lagrange multipliers vis), and \{s) can be identified as the x 
and y components of the internal force: Hx^v and ny = A. Requiring (3.13) to vanish for all x 
yields, after integration by parts: 



{-Hw {y{D) -y)-rjw + n) x 



D 



rD 



[nwy' + v'^ xds + 



UeX 



v'^xds = Q 



The fact that we have x(0) = 0, but arbitrary x{D) and x{L) implies: 

fc'e(L)=0, VeiD) - VilD) = -rjW , i^e = 0> ^'i = -IJ'Wy' 
Requiring (3.13) to vanish for all y similarly implies: 

Ae(L) = 0, \e{D)-\i{D) = -iiwx{D), A'g = 0, >!i = ^iwx' 



(3.22) 



(3.23) 



(3.24) 



where we see that Laplace pressure generates an outward normal force fiw{y' , —x') that causes 
the internal force (i^j, Aj) to vary. In addition we see that the internal force vanishes at the s = L 
extremity and that it experiences a discontinuity at s = D. Using (3.4) and (3.16) we find that: 



cos /3 
— sin/3 



(3.25) 



that is the external force applied on the beam at s = D is along the meniscus, as used in Neukirch 
etal. (2007); Antkowiak e/ a/. (2011). 

Once 7, A7, YI, and the volume V are set, the equilibrium configuration is found by 
solving the nonlinear boundary value problem for s G (0; D), given by Eqs. (3.5), (3.20), (3.21), 
(3.23), and (3.24), with left boundary conditions (3.7) and right boundary conditions 6'{D) = 0, 
= — jw cos /3, and \i{D) = jw sin (3. The presence of unknown parameters (5, r, and D 
is balanced by additional conditions (3.3), (3.4), and (3.18). 

We remark that, as in the case of a rigid substrate (Section 2), the sum of the distributed 

Laplace pressure /i = —iiwx'ds applied along the y axis on the beam (Eq. (3.24)) is balanced 
by the y component of the meniscus force at s = D: /2 = sin (5, that is /i + /2 = 0. 



Equilibrium solutions 

We now solve the boundary value problem for different values of the parameters, e.g. A, 7, 

A7. We note 6y the wetting angle, defined by A7 + 7 cos Oy = 0, and we use By instead of A7 
as parameter. We start with non-dimensionalizing the equilibrium equations. The configuration 
of the beam in the region s G [D] L] is trivial: the beam is straight and there is no stress 9'{s) = 
M{s) = and Ve{s) = Ae(s) = 0. The value of the length L is therefore of no importance, it can 
be anything as long d& L> D. Consequently we use ^fA as unit length, and YI/A as unit force. 
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The equilibrium equations for the dimensionless quantities (with over-tildes) read: 

x'{s) = cos9, y'{s) = sine (3.26a) 

61" (s) = hx sin 0-hy cos 9 (3.26b) 

n^(s) = -(7/f) sin 6*, hy{s) = (7/f) cos9 (3.26c) 

where 7 = A^fw/{YI) is a dimensionless quantity measuring the strength of surface tension. 
Volume conservation (3.3) now reads: 



1 



^2 



^- ^sin2^ ]+2xDyD-2 



D 



y cos 9 ds 



(3.27) 



The boundary conditions at s = are x{0) = 0, y(0) = 0, and ^(0) = 0. The boundary conditions 
at s = -D are: 

x{D) = fsm/3, 9'{D) = 0, n^(L») = -7 cos /3 , fiy(D) =^ sin (3 , 0{D) + P = 9y (3.28) 

For each value of the fixed parameters 7 and 6y, we numerically solve this boundary value 
problem with a shooting method where the six unknovms ^'(0), nx{0), ny{0), (3, D, and f 
are balanced by the five boundary conditions (3.28) and the constraint (3.27). Results for the 
inclination of the beam at s = D are plotted in Fig. 3-left. We remark that Eqs. (3.26b) and 




e{D)if(eY) 




.0 0.5 1.0 1.5 2.0 2.5 3.0 



Aywl (77) 



Figure 3. Inclination 9{D) of a beam deformed by a liquid drop, in the absence of gravity. Left: 6{D) as a function 
of the non-dimensionaUzed surface tension 7, for several value of the wetting angle 6y. Right: 9{D) scaled with 
fiOv), see formula (3.30). Near the origin, curves for different 6y collapse on a straight line of slope unity. 



(3.26c) can be simplified to: 



T 2 



(3.29) 



and boundary conditions (3.28) imply (0) = 0. 

We also look for analytical solutions when the surface tension is small 7^1, that is when 
the elasto-capillary length is large: L^c = \/YTJ(yw) S> \fA. We develop unknowns in power 

of 7: i(s) = xq{s) + 7x1(5) + ... , y{s) = + 7yi(s) + • • • . 0{s) = + 76)1(5) + ..., i^* = 
Do + jDi + ..., etc. We find xo{s) = s, 9i{s) = s{3D'^ - s^)/{6r) and yi{s) = s'^{6D'^ - 
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^)/(24f), /3o = Oy, Dq = fo sin 6'y, and 1/fo = \JQy — sin^y cos^y. This yields: 

d t'y — sm C^y COS Vy 

57 sin'^^y ^/~2\ 



(3.30) 
(3.31) 



4. Liquid drop on a flexible beam, in the presence of gravity 



s=L 




Figure 4. A flexible beam bent by a liquid drop, in the presence of gravity. Due to symmetry, we only consider 

positive s's. 



We now consider the situation where gravity is no longer neglected. The beam has mass per 
length r and the weight of the beam introduces the term: 



vL 



Tgy{s)ds 



(4.1) 



in the potential energy (3.10). Yet the main difference is that the liquid-air interface is no longer 
a circular-arc, and its shape yiy{x) has to be solved for. The liquid-air interface energy is now: 



rx{D) 



(4.2) 



In addition the weight of the liquid (of density p) has to be accounted for, yielding the term: 



rxiD) 



pgw 



rviv 



y dy dx 



Finally, the constraint of constant volume now reads: 



■x{D) 




yiv 



dy dx 



(4.3) 



(4.4) 
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As done previously, the beam geometry is parametrized with the arc-length s: x{s), y{s), 
9{s). As for the liquid-air interface, we introduce the relative height h = yi^ — y{D) and we 
parametrize it with the same variable s: h{s) = h{x{s)), see Fig. 4. The arc-length s along the 
beam is thus the unique independent variable. Dropping out constant terms in the energy, the 
Lagrangian of the system reads: 



ji:{x,y,e,h,D) = ^YI 



L 



Tgyds + wD A'y + 



+ pgw 

rD 



D 





2 ..2 



D 



Vx'^ + h'^ ds 



{h + y{D)) —y x'ds — fiw 



D 



{h + y{D) - y)x'ds 



+ 



+ 







t'j(s) \_x' — cos9\ ds + t'e(s) \x' — cos9\ ds 



D 
rL 



\i{s) [y' - sin 61] ds -|- Ae(s) [y' - sin 6*] ds 

JD 



(4.5) 



Boundary conditions are a;(0) = y(0) = 0(0) = 0, h'{0) = 0, and h{D) = 0. As in Section 3, we 
note X = {x, y, 9, h, D) and we look for Xe such that the first variation C'{Xp^{X) vanishes 
when X is changed from X^. to Xe + eX. We compute the first variation of C with regard to the 
variable X = (x, y, 9, h, D): 



C'{X,){X) = YI 



L _ rL 

9' 9'ds + 



Tgyds + wD A'y 



x'x' + h'h' 

V^^Tw 



ds + D^w^x'{Df + h'{DY 



rD 



+ pgw 



hhx' + hy{D)x' + hy'{D)Dx' + hy{D)x' + y{D)y{D)x'+ 

x' > ds 



Dy{D)y'{D)x' - yyx' + - {h + y{D)f - y^ 



+ 



jiw 

'D 


•D 



D 



+ 



[{h + y{D) + Dy'{D) - y) x' + {h + y{D) - y) x'] ds 



Vi{s) [x' -I- 6' sin 6*] ds -|- Ve{s) [x' -|- 6* sin 6*] ds 

JD 

- - 
\i (s) [y' - 9 cos 9] ds + Ae (s) [y' - 9 cos 9] ds 

JD 



(4.6) 



We now require the first variation to vanish for all x, y, 0, h, and D. Collecting terms involving 
h and h' yields, after integration by parts: 
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The first term is 'ywh' {D)h{D) / -^f x' {DY + h'{D)^. Boundary conditions require that h + eh 
vanishes at s = D + eD; this yields h{D) = —Dh' (D). Consequently this first term effectively 
goes into the equation for D, see (4.10). The second term of (4.7) has hence to vanish for all 
h{s), which imphes that the hquid-air interface h{s) obeys the equation: 

[pgwih + y{D))- fxw]x' = -iw(^-^=^==^ (4.8) 

Integrating this equation from s = to s = D yields: 

pgwA — iJLwx{D) = sin ip{D) (4.9) 

where A = {h + y{D))x'ds is the area between the hquid-air interface and the horizontal 
axis, and where (p is the angle the interface makes with the horizontal. Evaluating (4.8) at s = 
reveals that the Lagrange multiplier is the hydrostatic pressure at the origin. 
Requiring (4.6) to vanish for all D yields: 

wAj + w-i^x'{DY + h'{Df - wj^=^^£=== + y'{D) (pgwA - fxwx{D)) = 

(4.10) 

Using x'{p)l^x'{pY + h!{Bf = cos (p{D), x'{D) = cos e{D), y'{D) = sin e{D), and (4.9) 
we arrive at: 

A7 + 7cos[0(D) + ^] = O (4.11) 

where (3 = —(p{D). This is Young-Dupre relation for the wetting angle 9{D) + p. We collect 
terms involving x and x' in (4.6) and we integrate by parts to obtain: 



D 

+ T.Pgw [h + y{D)Y - ^pgwy^ + i^w{h + y{D) -y)\x 



w-/x' , 1 _ , ^JT^^^2 1 2 



J 



+ pgw{h + y{D))h' - pgwyy' - pw{h' - y') + f -j xds 



+ 



u'^x<ls = Q (4.12) 

D 



I VVx'2 + /l'2 

r - 1 ^ 

u^x 

L \d 

The fact that we have x(0) = 0, but arbitrary x{D) and x{L) implies: 

ye{L) = (4.13) 

v'e{s) = Q (4.14) 

Ve{D) - Ui{D) = w-^x{D)/^/x'{Df + h'{Df = w-^ cos ip{D) = w-f cos /3 (4.15) 

^'i = - [iJ-w - pgwy) y' - [pgw{h + y{D)) - p.w] h' - (^-j=^==^ (4.16) 

Considering the identity x' (^x'/Vx''^ + h'"^^ + h' ^h' / y/x^^~+W^ = and (4.8), equation 
(4.16) reduces to: 

u'i = - (pw - pgwy) y' (4. 17) 
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Requiring (4.6) to vanisli for all y similarly yields, after the use of (4.9): 

Ae(L) = 0, Xe{D)-Xi{D) = -wj sin 13 (4.18) 

K = T9 , K = rg + ii^w - pwgy)x' (4.19) 

Finally, requiring (4.6) to vanish for all 6 yields the same equation as before, see Eqs. (3.20) and 
(3.21). 

From (4.15) and (4.18), we see that the internal force experience the same discontinuity as 
in the case without gravity (3.25): here also the external force appUed on the beam at s = D is 
along the meniscus. 

Equilibrium solutions 

We now solve the boundary value problem for different values of the parameters, e.g. A, 7, 
A7, r. We note 9y the wetting angle, defined by A7 + 7 cos 9y = 0, and we use By instead 
of A7 as parameter. We start with non-dimensionalizing the equilibrium equations. As the 
configuration of the beam in the region s € {D; L] is no longer trivial, we use L as unit length, 
EI/L'^ as unit force, and EI/L as unit moment. For the beam, the equilibrium equations for the 
dimensionless quantities (with over-tildes) read: 

x'{s) = cose, y'{s) = sine (4.20a) 
e" {s)=hx sin e-riy cos e (4.20b) 

n^(s) = -P sine, fi'y{s) = f + P cos 6* (4.20c) 

where P is the dimensionless hydrostatic pressure P = {L/Lec)'^ [fi — {L/Lcf y] for the region 
s G [0; D) and P = for s > JD. We have introduced the dimensionless pressure jl = nL/'~^, the 
capillary length Lc = \/"i/{pg), and the elasto-capillary length Lgc = \/YI/{'^w) (Bico et al, 
2004). 

The equations for the liquid-air interface (4.8) can be rewritten using (i) the angle if the 
interface does with the horizontal, and {ii) the arc-length a along this interface: 



^'{a) = {L/L,Y h + y{D) 



jl (4.21a) 



h'{a)=sin<f, = cos (/? (4.21b) 

where a = a jL. The liquid-air interface has total contour length I. 

As soon as values for the fixed parameter Lgc/ -Z^, -Z^c/ and By are given, the boundary 
value problem can be solved with a shooting procedure where ^'(0), na;(0), ny(0), /i, /i(0), I), 
t are seven unknowns. Integration of (4.20) is performed for s G [0; ^) and integration of (4.21) 
is performed for a G [0; £]. At s = L> a jump in the force vector is introduced according to (4.15) 
and (4.18) with /3 = —<^(€) and = u, Uy = X. Then integration of (4.20) is performed for 
s G {D; 1). Seven boundary equations have then to be fulfilled: 

^'(1) = , n,(l) = , hy{l) = , h{£) = , i{£) = x{D) , e{D) +13 = By (4.22) 
together with the volume condition: 

ri rD r- 



A: 



h cos ip da 





y - y{D) 



cos9ds (4.23) 
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A solution for Lec/L = 0.175, Lc/L = 0.982, f = 2,A = 0.039, and 9y = 110 deg. is shown 
in Fig. 5. The seven unknowns of the shooting procedure are found to be 6'{0) = 1.887, nx.(0) = 

y 



0.3 








0.2 








0.1 








0.0 








-0.1 


0.0 0.2 0.4 


0.6 0.8 1.0 


Figure 5. A flexible beam bent by 


a liquid drop, in the presence of gravity. 



X 



9.15, %(0) = -3.32, /i = 5.24, h{0) = 0.219, D = 0.201, and £ = 0.333. 



5. Pinning, receding, and advancing of the contact line 

We now turn to the case of contact Une pinning and we show, following Alberti and DeSimone 
(2011), that pinning, receding, and advancing of the contact line can be treated in a variational 
approach where irreversibility conditions are introduced. We consider the drop of Fig. 1, lying 
at equilibrium on a rigid substrate with its contact angle /3 equal to the Young-Dupre angle 9y, 
defined by A7 + 7 cos Oy = 0, /3(t = 0) = 9y- At time t = evaporation starts to take place 
and we study the subsequent behavior of the drop. In the absence of contact line pinning, the 
contact angle (3 will stay at P{t) = Oy, and the wetted length D will decrease in order to keep on 
fulfilling the volume constraint. In the presence of contact line pinning, the length D will first 
stay fixed and the contact angle f3 will decrease: /3{t >0) <6y to fulfill the volume constraint. 
Eventually, as (3{t) reaches a receding threshold, the contact line will start to move: D = D{t). 

As in section 2 we start with the interfaces energies, Eq. 2.1, and we utilize the constraint 
D = r sin (3 to eliminate the variable r: per unit w the energy is 27D/3/ sin /3 + 2DA7. As D 
decreases from D{t = 0) receding of the contact line is associated with an energy dissipation 
k>0 per unit area (Joanny and de Gennes, 1984; Alberti and DeSimone, 2005), which we 
introduce in the energy: 

E{D,l3) = 2jD^ + 2DA-f + k{D{0)- D) (5.1) 
sm p 

We minimize this energy under (i) the constraint of fixed volume V/w = r'^{P — sin /3 cos /3) = 
{D/ sin/3)^(/3 — sin/3 cos/3), and (ii) the irreversibility condition D{t^) < D{t~), which we 
note D < . We therefore introduce the Lagrangian: 

£ = E{D, /3) - X{D- -D)- f,D' - ^) (5.2) 

V sm^ p sm /3 / 

where /j, is the Lagrange multiplier associated with the volume equality constraint and where 
A is the Lagrange multiplier associated with the inequality constraint D < D^. The necessary 
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conditions for having a minimum are: 
dC 

— =O^^D = ^smP (5.3) 
= ^ 27cos/3 + 2A7-A; + A = (5.4) 

oD 

K.T. ^ A > , (L>- - L>) > , and \{D- - £>) = (5.5) 



where the last line lists the classical Kuhn-Tucker conditions (Luenberger, 1973) arising in 
case of inequality constraints. These three conditions express the fact that either the pinning 
force A is zero and sliding occurs D<D^, or the pinning force is strictly positive A>0 
preventing the contact line to move: D = D~ . We now introduce an angle 9* such that 
2(A7 + 7 cos ^?*) = fc. Positivity of the dissipation k implies that 9*<9y- Equation (5.4) 
becomes 27(cos (3 — cos ^*) + A = and ehminating A in (5.5), we finally obtain: 

cos 6''' - cos /3 > (5.6a) 

D- - D>0 (5.6b) 

(cos 9* - cos P){D- -D) = (5.6c) 

which means that either D is fixed to D~ and /3 > 9* (contact line pinning) or 13 = 9* and D 
decreases D < D (contact line sliding, here receding), as was used in Rivetti and Neukirch 
(2012). 

Note that for simplicity we have only presented equations for receding of the contact line, 
but the present treatment can be done for the general case where advancing and receding can 
both occur, see FedeU et al. (201 1). 

Illustration 

We here illustrate the present theory on an imaginary experiment where one deposits a drop 
on an elastic strip (see Fig. 4) and wait for evaporation to take place (Py et al, 2007). After 
deposition the drop contact angle 9{D) + (3 takes some intermediate value between receding 
(9*) and advancing values. We fix parameters L^c = 0.2 L, Lc = 0.8 L, and f = 1.4 and we first 
solve equations (4.20), (4.21), (4.22), and (4.23) for several values of the receding angle 9* in 
the sliding hypothesis 9{D) + 13 = 9* . We then solve the equations for several values of -D(O) in 
the pinning hypothesis D = D{0). Results are shown in Fig. 6. 

In a typical experiment, starting at ^ = A/L'^ = 0.14 with D{0)/L = 0.32, evaporation first 
results in the decrease of the contact angle, following the curve Cdi in Fig. 7. As the contact 
angle 9(D) + /? reaches 9* (with say 9* = 2) the system switches branch at point Pi and follows 

the constant contact angle curve Cq down to ^ = 0. If now one starts at ^ = 0. 14 with D{0)/L = 
0.4, during evaporation the system follows curve Cd2 down to point P2. The constant contact 
angle curve at P2 being unstable the system jumps to a configuration with same volume, not 
discussed here. 



6. Conclusion 

In conclusion we have recalled that the classical Young-Dupre relation for the contact angle of 
a drop lying on a rigid substrate can be derived from a variational approach where the concept 
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X(L) / L 




Figure 6. Bifurcation diagram for the system of Fig. 4 witii L^c = 0.2 L, = 0.8 L, and f = 1.4. Tiie first set of 
eleven (red) curves correspond to equilibrium in the sliding hypothesis 6{D) + P = 8*, with 6* — 1.2, 1.3, . . ., 2.2. 
These curves all have a limit point for the variable A — A/L^. Dashed parts of the curves correspond to unstable 
equilibria. The second set of ten (black) curves correspond to equilibrium in the pinning hypothesis D — D{0), with 
D{0) /L = 0.2, 0.3, 0.32, 0.34, 0.35, 0.36, 0.37, 0.38, 0.4, 0.45. 



of force is not invoked. The variational approach has then been extended to the case where the 
substrate is a flexible beam and we have shown that (i) Young-Dupre relation still holds, and that 
(a) the external force applied on the elastic beam at the triple point is tangential to the liquid- 
vapor interface. We then extended the approach to the case where gravity is included and found 
that these two results continue to hold. Finally we have illustrated our model with the study of 
the evaporation of a drop deposited on a flexible beam and we have shown that, depending on 
the initial spreading of the drop on the beam, evaporation leads to a flat or a folded system. 

The present result showing that in the case a flexural deformations the external force on the 
elastic beam is along the meniscus is in contradiction to what is found in Marchand et al. (2012) 
in the case of extensional deformation. We are now working on the extension of the present 
variational approach to the case of extensional deformations and hope to resolve this apparent 
disagreement. 
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Appendix 1: Planar elastica 



The equilibrium equations for the system in Fig. A.8 are (Love, 1944; Antman, 2004): 

x'{s) = cos9, y'{s) = smO (A.la) 



Yie'{s) = M, M'{s) = sin 9 - Ny cos 



Kis) 



(A.lb) 
(A.lc) 



where s is the arc-length of the beam, M the internal bending moment, N = {N^, Ny) the 
internal force, (x, y) the current position of the central line, and 9 the angle between the tangent 
of the centre line and the horizontal axis. The bending moment is hnearly related to the curvature 
9'{s) through the bending rigidity YI, where Y is Young's modulus and I is the second moment 
of area of the beam cross section. In the case of a rectangular cross section of thickness h and 
width w, I = h^w/12 when bending occurs in the plane of the thickness h. The beam is clamped 
at s = and a vertical force (0, — F) is applied at the s = L extremity. We also consider the self- 
weight of the beam p'^^* = (0, —rg), where r is the mass per unit arc-length of the beam. The 
left boundary conditions are x(0) = 0, y(0) = 0, 9{Q) = 0, and the right boundary conditions are 
N^{L) = 0, Ny{L) = -F, and M(L) = 0. 

The equiUbrium equations can be recovered by considering the energy: 



E=^YI 
2 



[9fds + rg 



yds + Fy{L) 



and the Lagrangian: 



Cix,y,9) = E + 



{s) \x' — cos 9'\ ds + \i{s) [y' — sin ^] ds 
Jo 



(A.2) 



(A.3) 



subjected to the left boundary conditions. The conditions for the vanishing of the first variation 
of the Lagrangian will yield the equilibrium equations (A.l) together with the right boundary 
conditions (Audoly and Pomeau, 2010). 
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Figure 7. Evaporation experiment with Lee = 0.21/, I/c = 0.8I/, and f = 1.4. The curve Ce corresponds to 
equilibrium in the sliding hypothesis 9{D) + f] = 6* with 6* — 2; the region 6{D) + /3 < 2 is shown shaded. The 
curves Cdi and Cd2 correspond to equilibrium in the pinning hypothesis D = D{0), with D{0)/L = 0.32, and 0.4 
respectively. If evaporation starts on Cdi at A/ L"^ = 0.14 the system reaches point Pi and then follows the upper 
part of curve Ce down io A — Q. If evaporation starts on Cd2 at A/L^ = 0.14 the system reaches point P2. As 
configurations on Ce are unstable in this region, the system jumps to a configuration with the same value of A, to be 
discussed elsewhere. 
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Figure A.8. Cantilever beam sagging under the combined actions of its own weight r g and of a localized shear force 

F at its left end. 



